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Abstract. Helicases are molecular motors which unwind double-stranded nucleic acids (dsNA) in cells. 
Many helicases move with directional bias on single-stranded (ss) nucleic acids, and couple their directional 
translocation to strand separation. A model of the coupling between translocation and unwinding uses an 
, interaction potential to represent passive and active helicase mechanisms. A passive helicase must wait 

O |i for thermal fluctuations to open dsNA base pairs before it can advance and inhibit NA closing. An active 

helicase directly destabilizes dsNA base pairs, accelerating the opening rate. Here we extend this model to 
include helicase unbinding from the nucleic-acid strand. The helicase processivity depends on the form of 
the interaction potential. A passive helicase has a mean attachment time which does not change between 
ss translocation and ds unwinding, while an active helicase in general shows a decrease in attachment time 
xjy i during unwinding relative to ss translocation. In addition, we describe how helicase unwinding velocity and 

O ' processivity vary if the base-pair binding free energy is changed. 
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PACS numbers: 82.39.-k,87.10.+e,05.40.-a,82.20.-w,87.15.Aa,87.15.Rn 

1. Introduction 



00 , Helicases are motor proteins which separate the two strands of helical double-stranded nucleic acids (NA). 

£T) ' Both DNA, RNA, and DNA-RNA hybrid helicases are found in cells. Strand separation requires breaking 

the base-pairing interactions between the two strands and therefore requires energy input. Unwinding is 
fueled by NTP hydrolysis, typically of ATP. Helicases play a role in nearly every cellular process which 
involves NA, including DNA replication and repair, recombination, transcription, translation, and RNA 



processing^]. Aberrant functioning of helicases is associated with genome instability (the accumulation of 
*"^» , damage and errors in the genome), premature aging, and cancer 

q ' The essential common feature of all helicases is their ability to move along NA strands and couple motion 

to strand separation. (We use the terms opening, strand separation, and unwinding interchangeably.) For 
this reason, helicases are also NA translocases and share some features with other proteins which move on 
! NA strands 011]. 

The unwinding velocity and processivity of helicases are important for helicase function. The unwinding 
velocities of helicases range from tens to thousands of base pairs per second|T]. Helicase processivity is also 
variable. Processivity is most often defined as the average number of base pairs unwound per helicase binding 
event. Measured values of helicase processivity range from tens to tens of thousands of base pairs^Q. The 
processivity can be significantly altered by accessory proteins El Uii polymerasesjH] and multiple copies 
of the same helicase|??l 1101 HTj . Attachment of an enzyme either to a surface or another protein tends to 
increase the processivity ^2]- I n single-molecule experiments on RecBCD helicase, a lower processivity was 
measured for a helicase that was free in solution^] than for a helicase attached to the surface [TJ. The 
velocity and processivity of helicase proteins arc thought to be related to their biological roles. For example, 
replicative helicases are responsible for unwinding all cellular DNA during DNA replication. It has been 
proposed that replicative helicases should therefore have a high velocity and processivity, so that a small 
number of helicases can efficiently function in DNA copying. By contrast, a helicase which functions in DNA 
damage repair may only need to unwind a small region of DNA near a damage site; therefore its velocity 
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Figure 1. Sketch of helicase on a nucleic acid strand. The helicase moves forward — toward the junction 
where the dsNA strand opens — at rate fc+ and backward at rate k~ . The NA opens at rate a and closes at 
rate f). 



and processivity may be low. Understanding the physical basis of the velocity and processivity will help 
illuminate how these proteins may be optimized for different cellular roles. 

Different definitions of processivity have been given in the literature. They include (i) the average 
time a motor stays attached to its track, (n) the average number of steps a motor travels before detaching, 
(m) the number of ATP molecules hydrolyzed before the motor falls off the track, and (iv) the probability 
that a motor takes one more step (as opposed to unbinding before the next step). For helicases and other 
motors which interact with obstacles, this picture is complicated. We may consider both the translocation 
processivity, which is the average number of forward steps taken by the helicase taken during one binding 
event, and the unwinding processivity, which is the average number of NA base pairs (bp) unwound during 
a helicase binding event. The translocation and unwinding processivity are not equivalent, and their values 
vary with initial conditions. Unwinding processivity is usually measured for helicases; however, the measured 
value can depend on where the motor binds relative to where it begins unwinding the NA. 

Because helicase proteins are involved in a wide range of cellular processes, many different types of 
helicase protein are found in cells. Even the relatively simple bacterium E. coli has at least 11 different 
types of helicase. Helicases are structurally diverse, and not all helicases share a common mechanism. 
However, many helicases share the ability to translocate directionally on single-stranded (ss) NA, especially 
members of superfamilies I and 110] . Motion on ssNA is analogous to the motion of a classical motor protein 
which moves with directional bias on an infinite, one-dimensional lattice (figure 1). Single-stranded NA is 
polar, with one end labeled the 3' end and the other the 5' end. Helicases which directionally translocate on 
ssNA are referred to as 3' — > 5' helicases or 5' — > 3' helicases, depending on which direction they move on 
the single strand. When the helicase is near the ss-double strand (ds) junction, the helicase can move the 
junction forward, creating additional ssNA "track" as it moves. 

A natural question in the study of helicases is how a protein may efficiently couple translocation to 
unwinding. In the biochemical literature on helicases, this coupling is classified as passive or active [T51 HI IHl . 
A passive helicase waits for a thermal fluctuation that opens part of the dsNA, and then moves forward, 
binds to the newly available ssNA, and prevents the NA from closing. An active helicase directly destabilizes 
the dsNA, presumably by changing the free energy of the ds state. 

Many different types of helicase have been studied experimentally in some detailjU Ej- In recent 
years, single-molecule experiments have been performed for several helicases, directly measuring helicase 
unwinding rates. The helicase RecBCD, which is also an exonuclease (it degrades one of the unwound 
strands as it moves) has been studied by several groups [T71 1131 IT51 ITU ITU] . Interpretation of the RecBCD 
data is complicated by the fact that RecBCD contains two helicases of opposite polarity, which move on 
both strands of the DNA|20l 1211 122j . Bulk kinetic experiments on RecBCD have tried to determine the 
kinetic step size by fitting the data to a model of helicase unwinding [2*31 l2~H l2~5l 1215] : to date physical steps of 
this helicase have not been directly observed. Single-strand translocation and double-strand unwinding were 
observed for Rep helicase using single-molecule FRET [771 I28| . The motion of the RNA helicase DbpA has 
been observed with AFM[2*§]. Single- molecule experiments on UvrD helicase have observed unwinding, and 
evidence of strand-switching 30 . Finally, single-molecule FRET has been used to study the T4 replisome, 
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the replication machinery which includes the T4 helicase 31 . 

Although extensive biochemical and structural studies of helicases have been performed, few descriptions 
of the physics of helicase unwinding exist. Previous work includes a "flashing field" model specific to 
hexameric ring helicases 32 j, a description of a helicase as a biased random walk, which considered how 
the density of histones affects the random walk|33j , and a model where helicase motion is represented as a 
propagating front |34l I35j . Molecular dynamics simulations have addressed properties of Per A helicase 36 . 
Recent work includes a proposal that HCV helicase functions by a ratchet mechanism [5?J, a theory for 
two coupled motor proteins, proposed to describe RecBCD helicase |3*5|. and a detailed description of the 
mechanochemistry of T7 helicaseO- A physical description of helicase unwinding of NA has been proposed 
which contains both active and passive opening as different cases in a general framework |4Hl l%T) . Here we 
extend the description of helicase unwinding to allow calculation of helicase processivity. 

We represent a helicase by a particle which moves with directional bias on a one-dimensional lattice 
(corresponding to the ssNA). The ss-ds junction represents an obstacle on the lattice which blocks the helicase 
motion while also moving due to thermal fluctuations. Our model describes different forms of the interaction 
between the particle and the mobile obstacle. This scenario — a moving particle which interacts with an 
obstacle — also exists in other biological systems. Peskin, Odell, and Oster introduced the "polymerization 
ratchet" to describe how a growing biological polymer can exert a force against a fluctuating obstacle 02 • 
They argued that the rate of polymerization is limited by the time required for the obstacle to diffuse one 
monomer size. In the language of our model, described below, this scenario corresponds to a hard-wall 
interaction potential between the growing tip and the obstacle. As we show here, other forms of interaction 
can show significantly different velocity and processivity. Another example of the interaction of a motor 
with a second degree of freedom is the kinesin- family motor protein MCAK, which induces microtubule 
depolymerization if it interacts with the microtubule endgnHHESl. This process is relevant in the generation 
of force during the separation of chromosomes by the mitotic spindle. Finally, the interaction of two motor 
proteins on a filament is another example for which our approach is relevant. 

In this paper, we first review a simple description which captures both active and passive unwinding by 
helicases0U|E]- This framework considers a single helicase which does not unbind from the NA strand, and 
permits us to calculate the unwinding velocity for a given interaction potential (section [3J|. The unwinding 
velocity for a passive helicase with hard- wall interaction potential is typically significantly slower than the ss 
translocation rate of the motor far from the ss-ds junction. For a simple form of active opening, the velocity 
approaches the single-strand translocation rate of the motor. In other words, an optimized active helicase 
can unwind NA as fast as it translocates on ssNA. Comparable rates of ss translocation and unwinding is 
therefore a signature of active opening. 

We then extend this description to include a nonzero unbinding rate, and develop a simple model of 
unbinding that captures the key effects. We assume that the unbinding rate may be different for a helicase 
translocating on ssNA than for a helicase unwinding dsNA, because the interaction potential alters the 
unbinding rate when the helicase is near the ss-ds junction. If the interaction free energy is larger, than 
the helicase will unbind more quickly. In particular, we assume that the free energy difference between the 
bound state and the barrier to unbinding is decreased by the amount of the interaction potential. We assume 
that the unbinding rate is independent of the helicase biochemical state. The unbinding rate is determined 
only by the height of the energy barrier separating the bound from the unbound states. 

This model of unbinding allows us to calculate different measures of processivity (section |3J). We discuss 
the dependence of the processivity on the shape of the interaction potential. The average attachment time of 
a helicase is particularly simple for a passive helicase: the helicase unbinding rate during unwinding is equal 
to the unbinding rate during ss translocation (section[«OJ|. By contrast, the helicase attachment time is lower 
during unwinding than during ss translocation for an active helicase (section QJ. Decreased attachment time 
during unwinding is a signature of active unwinding. 

Finally, in our conclusion, we discuss the connection between our results and experiments (section |SJ). 
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2. Active and passive unwinding 



In our simplified discrete description of helicase unwinding, the position of the motor — the helicase — along 
its track is labeled by the integer n, and the position of the obstacle — the ss-ds junction — is denoted m 
(figure nj. The motor is assumed to move toward increasing n, and we expect n < m. The interaction 
between motor and obstacle is characterized by an interaction energy U(m — n) which depends only on the 
obstacle- motor separation. We assume U — ► for m 3> n and U — ► oo for m < n. 

One simple form of this potential is a hard-wall interaction of zero range: U = for m > n and U = oo 
for m < n. For a hard-wall potential, a motor which is near the obstacle (n = m — 1) can advance only if 
the obstacle undergoes a fluctuation which increases m. This situation corresponds to passive unwinding. 
For active unwinding, the interaction between the motor and obstacle has nonzero range: when the helicase 
is near the ss-ds junction, both the kinetics of NA opening and the helicase motion are altered by the 
interaction. We will discuss interaction potentials which represent active unwinding below. 

We denote the rates of forward and backward hopping of the motor on the lattice far from the obstacle 
by k + and k~ . The hopping rates for forward and backward motion of the junction (corresponding to the 
opening and closing of dsNA) in the absence of the helicase are denoted a and (3. The ratio of these rates is 
given by a/ (3 — e~ AG where AG is the free energy difference per base between dsNA and two complementary 
ssNA strands. (Note that we use units where fcgT =1.) This energy difference is positive when dsNA is 
thermodynamically stable. For simplicity, we write a similar expression for the ratio of forward and backward 
hopping rates of the helicase, k + /k~ = e AM . Here A/j, denotes the chemical free energy of ATP hydrolysis 
which drives helicase motion. This expression applies when each hydrolysis event is tightly coupled to a 
forward step on the NA. We use this simplification here to discuss principles which do not depend on the 
validity of this assumption. Note that an effective value of Afi can be derived in certain limits from a more 
detailed description of the mechanochemical coupling 41 . 

The interaction between the helicase and the ss-ds junction modifies the hopping rates. We express the 
ratios of these rates as 

Pj I 3 _-[uU-i)-u(j)] 



r e -l<J<J-V-UM, (la) 



k + k+ 

= _ e -[uu-i)-u(j)] , lh \ 
kj_ x k~ 

where otj, Pj, kj~ and k~ are the position-dependent rates when the helicase and the junction are separated 
by j = m — n bases. To fully specify the position-dependent rates we write 

k+=k + e-WU- 1 '*- u M, (2a) 

k - =k - e -(f-WU-V-uV)] t ( 2b ) 



P^Pe-WV- 1 )-"™, (2c) 

a j - 1 =ae-V- 1 X u V-V- u W. (2d) 

where the parameter / describes the energy barrier of the transitions. For a one-dimensional reaction, / 
corresponds to the fractional distance of the peak of the barrier between the two adjacent states and thus 
0</<l. 

For the simple case of a hard- wall potential (figure the helicase prevents NA closing when j = 1. 
Therefore Pi = because £7(0) is infinite. Infinite U(Q) also implies k^ = 0: the helicase must wait until a 
thermal fluctuation opens the NA before advancing. An interaction potential with nonzero range corresponds 
to enzymatically assisted opening. For simplicity, we use a linear potential with range N and slope U . The 
discretized potential is shown in figures[2p and|2t. The potential energy increases in N steps, each of energy 
Uq, before a hard wall is reached. The increase in energy due to this potential facilitates opening of the 
dsNA: U(j — 1) > U(j) implies pj/otj-i < P/a. In addition, the interaction energy slows helicase forward 
motion. We assume U — > for j — > oo, so that no interaction occurs for large separations. 

The probability P(j,l,t) that the helicase and junction are at separation j and midpoint position 
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Figure 2. Interaction energies between the helicase and ss-ds junction as a function of their distance m — n 
in number of bases, (a) Hard-wall potential, (b) Potential with a single step, (c) Potential with three steps. 



I = m + n satisfies the master equation 
dP(j I t) 

U df ' ' = - (ay + (3, + k+ + kJ)P(j, I) + aj-iP(j -1,1-1) 

+ j+1 P(j + 1,1 + 1) + k+ +1 P(j + 1,1-1) + kJ_ x P{j -1,1 + 1). (3) 

Since the rates depend only on j, we define the distribution Vj = P(j, I, t). After a relaxation time which 
depends on the rates a, (3, k + , and k~ , the distribution Vj relaxes to a stationary state which satisfies the 
recursion relation 

kj + a.; 

"Pj+i = -r V r ( 4 ) 

k j+l + Pj+l 

This steady state distribution can be used to calculate the mean velocity (bp s~ 1 ') |4(JlllT] 

r • k t ( 5 ) 

3 

This expression for v has a simple physical interpretation — the quantity in parentheses is the unwinding 
rate at separation j, which is multiplied by the probability Vj of finding the complex at separation j. The 
effective diffusion coefficient which characterizes velocity fluctuations is 

D =\Y,( k j+ a J + k J + ■ ( fi ) 

3 

We assume that the helicase remains bound to the NA for all time. Therefore the expressions for v and D 
are true steady-state values reached by a real system in the long-time limit. In the results below, unwinding 
velocities of a bound helicase are calculated under this steady-state assumption. 

3. Processive and unprocessive unwinding 

In section [5] we assumed that the motor never unbinds from the track, so the model represents a motor 
which is infinitely processive. Real motors have finite processivity — they unbind eventually. This effect can 
be incorporated in our description by introducing the rate of unbinding 7. We assume that 7 depends on 
the separation j between motor and obstacle, but has no other position or time dependence. For a repulsive 
potential such as the linear potential discussed above, when the motor and obstacle are close to each other 
the free energy of motor and junction is increased. This typically leads to an increased rate of unbinding. 
This effect can be described by writing 46 

7i = 7 e c/ «. (7) 

Here 7 > is the detachment rate for motion on a single strand far from the junction. 
The time evolution of the probability distribution is then described by 

= -(*:/ + kJ + ay + 0j + jj)P(j, I) + aj^PU -1,1-1) + f3 j+1 P(j + 1,1 + 1) 

+k+ +1 P(j + 1,1-1) + k-_ x P{ 3 -1,1 + I). (8) 
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Note that this expression differs from equation only in the additional term proportional to the unbinding 
rate jj. The other rates are unchanged as defined above, in equation (J2J. We assume that at time t — 0, 
the helicase and junction have a specific position determined by j = j a and / = l Q . Therefore, 

P(j,l,t = 0) = 6 jjo 6 Uo , (9) 

where Sjk = 1 for j = k and otherwise. 



3.1. Processivity and attachment time 

We characterize the processivity by the average attachment time (r) of the motor and by both the 
translocation processivity (Sn) , which is the average number of forward steps taken by the helicase during one 
binding event, and the unwinding processivity (8m), which is the average number of NA bp unwound during 
one helicase binding event. These measures of processivity depend on initial conditions. If the helicase binds 
near the ss-ds junction, then the translocation and unwinding processivities are similar. However, if the 
helicase binds to a ss region far from the junction, then the unwinding processivity can be much smaller 
than the translocation processivity. If the NA can close on average under these conditions, the corresponding 
unwinding processivity becomes negative. 
The mean attachment time is defined as 

/>OC 

(t)= dtty2 7j p(j,i,t). (io) 
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Similarly, the translocation and unwinding processivities can be defined by 

({m) _ w+M. (12) 

Here, the average changes in j and I during unwinding are 

p OG 

(8j)= / dt VPUM), (13) 

Jo 3,1 

p OO 

(si) = / dt Y,(i-io)-yoP(j,i,t)- (14) 

J ° 3,1 

The translocation processivity (8n) = (n — n ) is the mean number of steps the motor takes in one direction 
before detaching and the unwinding processivity (8m) = (m — m Q ) is the mean number of steps the obstacle 
moves per binding event. The above expressions can be simplified using Laplace transforms The Laplace- 
transform of P(j, I, t) is 

/•oc 

P jtl (s) = / dt e' st P(j,l,t). (15) 
Jo 

Summing equation © over j and I, we obtain 

£^^ = £-vwt). (i6) 

3,1 3,1 

We can thus simplify the expression for (r) in equation (|10J) by integrating by parts: 

dP(j,l,t) 



< t >=-e r dtt< - 

3,1 j ° 



dt 



p OO 

= V/ dtP(j,l,t). (17) 

3,1 J ° 

The boundary terms vanish because at the lower boundary t = 0, while for long times the probability 
approaches zero at all sites because the motor unbinds. 
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Defining the coefficients 

Qi,l = Pj,l(s = 0)= / dtP(j,l,t), (18) 
Jo 

the average attachment time can be expressed as 

(r)=X>.J" (19) 

3,1 

Similarly, the j and I processivities are given by 

/>oo 

(Sj)= / dt - Jo) 7jP(3,l,t), (20) 

J ° 3,1 

= 5> '- 3ohjQj,U (21) 

and 

(8l)=J2(l-loh 1 Q 1 A- (22) 

We therefore need to calculate Qj,i, the Laplace-transformed probability at a = 0, to obtain the average 
attachment time and the mean change in j and I during a single binding event. The translocation and 
unwinding processivities are obtained by equations 1)11(1 and (|12[1 . 

The coefficients Qjj are obtained by solving the Laplace-transform of equation JHJ for a = 0: 
' dP(j I) ~ 

dt e~ st g = -(At + fc. + aj . + /j. + 7 i )P ijI (s) + aj-iPj-i^-iCs) + ^ +1 P j+ i, i+ i(s) 

+kf +1 P j+ i,i-i(s) + kJ^Pj-iwia). (23) 
Note that the rates have no explicit time dependence. As above, integrating by parts leads to 

--st dP(j,l) r>f-j x\„-st|°° , „ f°° i, -st-. 



dt e~ st = P(j,l,t)e- St \^ +s J dt e~ st P(j,l,t), 

= -P(j,l,t = 0) + sP j , l (s), 

= -Sjj S tto +aP j , l {a). (24) 
Here we have used the initial condition, equation Thus the Laplace-transformed equations are 
~ S jjo S Uo + sp ],i{ s ) = ~( k t + k J + a 3 + + Jo)Pj,l{ s ) + otj-iPj-i,i-i(s) + f3 ]+1 P j+u+ i(s) 

+k+. 1 P j+1 , l - l (a) + kj^Pj-u+xis). (25) 
In general, these equations could be solved for arbitrary s. We only require the Laplace-transformed 
probability evaluated at s — 0. The coefficients Qj i satisfy the equation 

SjjJtta = + kJ + otj + Pj + 7j)Qj,j - aj-xQj-i,i-i - Pj+xQj+x,i+i - kf +1 Qj + x,i-i - kJ_ t Qj-x j+i-(26) 
These equations describe a discrete drift-plus-diffusion system in two dimensions with an inhomogeneous 
term. 

Equation i|26|) can be solved by a product ansatz 

Qj,i = RjTi (27) 
Using this ansatz in equation 12tj|) . and choosing the normalization Yli'Pi ~ 1 an( ^ ^j"fjPj = 1> the 
coefficients Rj satisfy 

" Sjj„ = -{k+ + kj + a, + 0j + lj )R j + (oj-! + kJ^R^x + (P j+ x + k+ +1 )R ]+1 . (28) 
The coefficients 1} satisfy the equation 

Su^ip + q + W-pT^x-ffi+x- (29) 
Here we have defined 

P = J2( a 3+k+)R 3 , (30) 

j 



3.2. Simplified expressions for processivity 

The average attachment time (t) and the quantities (6j) and (61) can be expressed in terms of the Rj and 
//: 

(r) =X)iJ,-, (32) 

j 

(Sj) = J2(J-Jo)1jRj, (33) 
I 

The solutions to equation (|28|l depend on the shape of the interaction potential. Some examples are 
discussed below. However, in general equation i|29|l can be solved formally. Equation (|29[1 is a second-order, 
linear, inhomogeneous difference equation. Since we have assumed translational invariance in I, the resulting 
solution is independent of l Q . Therefore we choose l = for convenience. The solutions have the form 
Ti = y l for I ^ 0. Using this ansatz, y obeys a quadratic equation 

y 2 -(l + a)y+{a-b) = 0, (35) 

where a — (l+p)/q and b — l/q. We denote by y+ and j/_ the positive and negative roots of equation i|35|l . 
For I > 0, Ti = A'y l _, while for I < the solution is X} = Ay l +1 since Ti must vanish for I — > ±oo. Requiring 
that To is the same for both expressions implies A — A'. Finally, equation l|29|l for I = leads to 

-6^ - 1 



A = b + a - y_ — 1 . (36) 

Formally, we can use equation l|34(l to determine (61): 

(6l) = ^2lT h (37) 

OO 

= Aj2{iy l --iy+), (38) 



= .4 



y- y-\ 



(39) 



[l-y-Y (l-y+) 

Note, however, that A and y± depend on the shape of the interaction potential via the coefficients Rj. Once 
these coefficients have been determined, the quantities (t), (6j), and (61), and therefore (6ri) and (6m), can 
easily be obtained. 

sectionVelocity and processivity for a hard-wall interaction potential 

3.3. Velocity 

For the hard-wall potential, the unwinding velocity is given by 

ak + - Pk~ 

v HW = - TT ^ T - (40) 

The velocity is positive whenever k + /k" > (3 /a, that is, the free energy change AG which drives NA 
closing must be smaller than the chemical potential A/i of ATP hydrolysis. The maximum vhw occurs for 
a unidirectional helicase (fc~ = 0). This upper bound is 

« - 1 (rnhfi) - (41 » 

where the approximation holds if k + <g; f3. Thus a passive helicase unwinds more slowly than it translocates 
on ssNA by a factor w a/ (3. This result has a simple interpretation: the base pair adjacent to the helicase 
has a probability a/j3 of being open. Thus, when the helicase attempts a forward hop it succeeds with 
probability a//3. 
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Figure 3. The translocation processivity (8n) and unwinding processivity (<5m) as a function of the initial 
separation j a between motor and junction. When j a = 1, the translocation and unwinding processivities are 
approximately equal. As the separation j increases, the translocation processivity (Sn) is approximately 
constant, while the unwinding processivity {Sm) decreases. The parameters are a = 10 5 s — 1 , j3 = 7 X 10 
s — 1 , fc + = 1 base s — 1 , k~ = 0.01 base s — x , 7 = 0.03 s -1 . Calculations were done numerically using a grid 
size of M = 100 in j. 



3.4- Processivity 

For a hard-wall potential, the average attachment time, translocation processivity, and unwinding 
processivity can be calculated analytically. We solve equation f2"%|l for Rj for the simple case j a = 1. 
For j > 1, the rates are independent of j and we have 

flj+i = (1 + d)Rj + (d - e)i?j_i, (42) 

where we have defined 

a + fc~ + 7 
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(43) 
(44) 



This equation has solutions of the form Rj = x±, where the x± are the positive and negative root of the 
quadratic equation x 2 — (1 + d)x + (e — d) — 0. The boundary condition that Rj vanishes for large j imposes 
Rj = Bx°_- From equation (gHJ at j = 1, we find (j3 + k + )(R 2 - dR x ) = 1 and 

B = [x-((3 + k+)(d-x-)}- 1 . (45) 

With these expressions, (r) = 7 _1 . This simple result follows since for the hard wall case jj = 7 is 
independent of the separation j. 
The average change in j is 

(8j) = "j ~~ (46) 
1 — X- 

To determine (SI), we note that for the hard wall case 

p = 7~ 1 (a + fc+a;_), (47) 
g = 7 _1 (fc" +)9ar_). (48) 

The quantity (SI) now follows from equation Ij39(l . In the hard- wall case, the translocation processivity is 
approximately the unwinding velocity times the attachment time: (Sm) w v(t), as we would expect for 
position- independent detachment. Note that the this result depends on initial conditions and is valid for 
./.. 1- 

When j = 1, the translocation and unwinding processivities are approximately equal. As the separation 
between the motor and obstacle increases, the unwinding processivity drops. A negative unwinding 
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Figure 4. Velocity and processivity as a function of the base-pairing energy AG for a hard wall-potential. 
We varied the opening rate of the ss-ds junction, assuming that the closing rate is constant. Therefore 
a = f)e~^ G . Left, unwinding rate vhw/ v ss as a function of AG. Right, translocation and unwinding 
processivity as a function of AG. Parameters are P = 7 x 10 5 s -1 , fc+ = 1 base s" 1 , k~ = 0.01 base s _1 , 
7 = 0.01 s — 1 , and j = 1. 



processivity reflects net backwards motion of the obstacle while the motor is bound. For a large separation 
between the motor and obstacle, they have no effect on each other. In figure yj the processivities are shown 
as a function of j a with other parameters fixed. 

The helicase adjacent to the junction can only move when the junction opens. Therefore, the velocity 
and the processivities depend strongly on the base pair binding free energy AG. The velocity and the 
measures of processivity are displayed in figure 0] as a function of AG. In this calculation, we assume that 
(i is constant and a — (3e~ AG . The value AG = 2 corresponds to the sequenced-average value for typical 
genomic DNA. For small AG, the processivity increases dramatically, because the binding energy driving 
NA closing is decreased. 



4. Velocity and processivity of an active helicase 

4-1- Velocity 

Active opening can be represented by an interaction potential with nonzero range between junction and 
helicase. For simplicity, we choose linear potentials characterized by a range of N steps and a step height U 
(figure|2J. The repulsive interaction between helicase and junction implies that for small separation j within 
the range of the potential the rates of junction opening and helicase backward hopping are increased. The 
unwinding velocity for such a potential, relative to the hard-wall case, is |4UI 14 1| 

v N _ c N (e- u ° - c) + (1 - c)e-f u °(e- NU ° - c N ) 

vhw ~ c N {e~ u ° - c) + (1 - c)e-Uo ( e -NU _ C W) ■ t 49 ) 

Note that the helicase cannot increase the unwinding rate beyond a certain limit. Assuming < / < 1, the 
unwinding rate has an upper bound vn < c~ 1 vhw- 

The unwinding velocity for one step (N — 1) is displayed in figure^ as a function of the step height 
U for different values of the parameter /. For small Z7o, the unwinding rate increases with increasing step 
height U because the presence of the step facilitates NA opening. For large step heights, the unwinding 
rate decreases, because the repulsive potential reduces the rate of helicase forward motion. The unwinding 
velocity for different values of N is shown in figure as a function of U . For increasing N, the opening 
rate becomes more sensitive to U Q and the maximum occurs at higher values. For large N, the maximum 
unwinding rate occurs for [/* ~ — lnc w AG. In this case fastest unwinding occurs when the slope of the 
potential matches the base-pairing energy of the NA. 
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Figure 5. Variation of unwinding velocity, average attachment time, and unwinding proccssivity as a 
function of step height for a one-step potential. The parameters are a = 10 5 b _1 , P = 7 X 10 5 s" 1 , k+ = 1 
base s — 1 , k~ = 0.01 base s — 1 , and, for the processivity calculation, 7 = 0.01 s" 1 , j = 1, and grid size 
M = 100. 
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Figure 6. Variation of velocity, average attachment time, and unwinding processivity with step height, for 
potentials with varying numbers of steps. The parameters are / = 0.05, a = 10 5 s" 1 , p = 7 x 10 5 s" 1 , 
fc + = 1 base s — 1 , fc~ = 0.01 base s — 1 , and, for the processivity calculation, 7 = 0.01 s — 1 , j = 1, and 
M = 100. 



4-2. Processivity for active unwinding 

Active unwinding results in a decrease of the mean attachment time, as compared to ss translocation or 
passive opening for which (r) = This is a consequence of the repulsive interaction potential which leads 
to increased unbinding rates for small j. The dependence of the average attachment time on the step height 
is displayed in figure [5Jd for N = 1. For a potential with one step, (r) decreases rapidly with step height to a 
plateau value w 0.27 -1 . This decrease is independent of /, since the unbinding rate does not depend on the 
barrier between states of different j. Note that the limits U and U — + 00 do not give the same result 
for the processivity because the initial conditions are different: t/ — > corresponds to a hard- wall potential 
with j = 1, while U — * 00 corresponds to a hard- wall potential but with j — 0. 

The decrease of (r) for increasing U is similar but even more pronounced for potentials with longer 
range, as shown in figure EJd. For N > 2, the average attachment time rapidly decreases to < 0.017 -1 . 
Decreased attachment time during unwinding (as compared to ss translocation) is a strong signature of 
active unwinding. 

Because the attachment time decreases with increasing step height, the translocation and unwinding 
processivities also tend to decrease. However, for a one-step potential (5n) and {5m) decrease slowly for 
increasing step height when / is small (figure Et). This occurs because the rapid increase in the velocity 
approximately cancels the decrease in attachment time, leading to a processivity that is relatively insensitive 
to step height. For larger numbers of steps, a decrease in (8m) occurs as U increases (figure Et) . These 
curves are approximately exponentially decreasing functions, with a decay constant determined by /. Small 
/ means that the processivity varies slowly in the biologically plausible range between — 20 fcsT. Note 
that above U a = 20, little or no increase in the velocity can be achieved by increasing the step size. 
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Figure 7. Unwinding velocity, attachment time, and unwinding processivity vs. free-energy per base pair 
for a one-step potential. We fix the closing and therefore a = f3e~^ G . The parameters are / = 0.05, 
/3 = 7 X 10 5 s _1 , fc + = 1 base s~ x , k~ = 0.01 base s — 1 , and, for the processivity curves, 7 = 0.01 s -1 , 
jo = 1, and M = 100. 



4-- 3. Dependence on base-pair binding energy 

The velocity and attachment time are sensitive to the mean free energy per base pair AG. We assume that 
the rate (3 is constant and a — f3e~ AG . For simplicity, we show results only for a potential with N = 1. The 
dependence of v and (r) on AG depends on the step height U . For a hard- wall potential, the velocity varies 
exponentially as a function of AG. As U increases, the velocity becomes less sensitive to AG (figure |7|i) . 
Thus for U = 5 the displayed curve decreases slowly with increasing AG. The attachment time (r) shows 
the opposite trend: for a hard-wall interaction potential, the attachment time is independent of AG. The 
attachment time is also weakly dependent on AG for small U , while for larger U Ql (r) varies more rapidly 
(figure [7|3). The processivity (5m) as a function of AG exhibits the same behavior independent of U (figure 
Et)j because the variation in velocity and attachment time as a function of step height approximately cancel. 
Therefore, measurements of velocity and attachment time have greater power to elucidate aspects of the 
potential. 

5. Discussion 

In this paper, we have extended a model for the unwinding of double-stranded nucleic acids by helicases to 
allow calculation of helicase processivity. An interaction potential describes how the helicase and NA ss-ds 
junction affect each other. We assume that the local unbinding rate of the helicase depends exponentially 
on the interaction potential. The mean attachment time, the translocation processivity, and the unwinding 
processivity depend on initial conditions: if the position where the helicase binds to the ssNA is varied, the 
processivity changes. Therefore, in different experimental situations different values of the processivity could 
be obtained. 

A hard-wall interaction potential describes passive unwinding. In this case, the unwinding velocity vhw 
is significantly slower than the single-strand translocation rate of the motor far from the ss-ds junction. The 
average attachment time for the hard- wall potential is (r) = 7 _1 since the unbinding rate 7 is independent 
of position. In other words, the mean attachment time of a passive helicase is the same whether the helicase 
translocates on ssNA or unwinds dsNA. For initial conditions where the helicase and junction are adjacent to 
each other, (5m) w vhw(t) — the translocation processivity is approximately the unwinding velocity times 
the attachment time. The translocation processivity is lower for a passive helicase during unwinding (relative 
to ss translocation); this decrease occurs solely because of the decrease in velocity of the helicase, while the 
average attachment time is unchanged. 

We represent active unwinding by a linearly increasing potential of finite range, characterized by the 
step height U~o and the number of steps N. A range of up to 10 bases is motivated by the typical size of a 
helicase monomer. Indeed, structures determined by X-ray diffraction suggest that a helicase can interact 
with 5-10 bases in the ds region of the ss-ds junction g3 ESI EH IS3 ED E21 EH • We find that for an active 
helicase the unwinding velocity can approach the single-strand translocation rate of the motor. In other 
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words, the helicase can unwind as fast as it translocates on ssNA. 

In contrast to the passive case, the attachment time of an active helicase is shorter during unwinding 
than during ss translocation. Decreased attachment time during unwinding is a general property of actively 
opening helicases in our model. Our representation of the unbinding rate assumes that the helicase is 
bound in a single potential well with a barrier separating the bound state from the unbound state. This 
model is an approximation, because the helicase exists in different biochemical states throughout the ATP 
hydrolysis cycle. Each state might experience a different free-energy barrier to unbinding. In principle, if 
the relative time spent in each of these states changes (for example, as the ATP concentration is varied), 
the unbinding rate might vary as well. However, experimental evidence suggests that our simplified model 
is a good approximation. Single-molecule measurements on RecBCDpS] and UvrD|3D] find that the average 
attachment times of these helicases are independent of ATP concentration. Therefore, the assumption that 
the potential well for a bound helicase is independent of hydrolysis state is consistent with experiments. 

In our model, an actively unwinding helicase shows a decrease in attachment time. Previous work has 
suggested that high helicase processivity may require two helicase-NA binding sites — one on the ssNA, and 
another on the dsNA[Tl 1531 155) . In this mechanism, a helicase may unbind from the ssNA but remain bound 
at the dsNA site. The opportunity for dsNA binding is absent when the helicase translocates on ssNA. 
This effect can be incorporated in our description by introducing an extra binding energy and consequently 
reduced detachment rate if the helicase is close to the junction. 

Both the velocity of unwinding and measures of processivity are sensitive to AG, the average free energy 
of NA base-pair opening. The value of AG can be controlled in single-molecule experiments where tension 
is applied to the ends of the NA, by the concentration of ss-binding proteins in the buffer solution, and by 
varying the base composition of the NA. 

Measurements of helicase velocity and processivity as a function of AG can provide information about 
the interaction between helicase and NA ss-ds junction. Because the rate of passive unwinding is determined 
by the opening probability of the base pair at the junction, the velocity and the processivity decrease rapidly 
as AG increases. For active opening, the unwinding velocity as a function of AG depends sensitively on the 
step height: while for larger values of U the velocity depends only weakly on AG, this dependence becomes 
strong for small Uq. The average attachment time (r) behaves differently. For a hard- wall interaction 
potential, the attachment time is independent of AG. For active unwinding, the attachment time depends 
strongly on AG for larger U Q and tends to decrease for increasing AG. The effects of AG on velocity and 
attachment time approximately cancel when computing the processivities, so that the behavior of (5m) as a 
function of AG is similar for different step heights. Galletto et al. found that the unwinding rate of DnaB 
helicase depends on the DNA GC content Varying the fraction of GC versus AT base pairs changes 
the average AG of the NA template. The dependence on GC content found in these experiments would be 
expected in our model for an interaction potential with a small (or zero) step height. 

For the helicase PcrA, crystal structures suggest that the protein binds both the ss and dsDNA and 
distorts the double helix [3] . If the protein residues which have been proposed to interact with dsDNA are 
mutated, the mutant proteins hydrolyze ATP in the presence of ssDNA at a rate similar to the wild-type 
protein, but unwind dsDNA 10-30 times more slowly than wild type [57] . In the language of our analysis, 
the mutations may alter the interaction potential between junction and helicase such that it resembles the 
passive (hard- wall) case. In our description, starting from a situation of active opening with optimally chosen 
step height and step number to a passive case with hard-wall potential typically leads to a decrease of the 
unwinding rate by a factor of cr 1 w 7. Further altering the potential to an attractive linear potential 40 
with a well depth of 2 /c#T decreases the unwinding rate by a factor of 35. Small changes to the interaction 
potential can thus cause the unwinding rate to vary by a large factor. 

In single-molecule experiments on UvrD helicase, Dessinges et al. observed the unwinding of DNA 
molecules by UvrD at different forces [HD]. In these experiments, the end-to-end extension of a tethered 
DNA molecule changes with time as a result of the transformation of dsDNA into ssDNA by the helicase. 
These single-molecule experiments found that unwinding can be induced by UvrD monomers. However, 
bulk experiments have suggested that helicase activity requires UvrD dimers|SSl|SHl|Hn], even though UvrD 
monomers have a large processivity (2400 bases) when translocating on ssDNA 60 . The reasons for this 
difference are not currently understood. The velocity of unwinding was shown by Dessinges et al. to depend 
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on ATP concentration|30|. In addition, events were observed where the DNA slowly closed at an ATP- 
dependent rate. It was suggested that these events occur when a UvrD molecule bound near the ss-ds 
junction switches from one DNA single strand to the other. Because the two single strands have opposite 
polarity, the helicase might move away from the junction after strand switching. In this interpretation, 
the re-zipping events provide information on ss translocation of UvrD. The re-zipping velocity would be 
expected to be greater than or equal to the ss translocation rate of UvrD, because the energetically favorable 
re-annealing of the two DNA strands behind the helicase might accelerate the protein's motion. 

We can relate three key measurements of UvrD motion by Dessinges et al. to the behavior of active 
unwinding discussed here. First, the experiments show that the re-zipping velocity and the unwinding 
velocity are comparable in certain situations: the re-zipping velocity is approximately 15% larger than the 
unwinding rate at 35 pN applied force. Second, the effective detachment rate is about 10 times larger 
during unwinding than during re-zipping (at 35 pN). Third, the measured unwinding velocity is only weakly 
dependent on AG. The velocity was measured for 3 pN and 35 pN applied force on the DNA. The first 
case corresponds to an increase of AG, while in the second case AG decreases |30|. If we interpret re-zipping 
events as approximating the ss translocation behavior of UvrD, all three of these experimental observations 
are consistent with our description of an active helicase. 

The comparison of these data to bulk measurements reveals agreement in the unbinding rate during 
ss translocation but not in the unwinding rate. The unbinding rate of UvrD monomers translocating on 
ssDNA measured in bulk by Fischer et al. is comparable to (2 times larger than) the unbinding rate 
measured by Dessinges et al. for UvrD monomers during re-zipping events 60 . Fischer et al. found that 
the UvrD unwinding rate was approximately 3 times slower than the ss translocation rate|6L)|. However, the 
bulk experiments of Fischer et al. observed unwinding only by UvrD dimers, whereas the single- molecule 
experiments of Dessinges et al. observed unwinding by UvrD monomers. Therefore the unwinding data may 
not be directly comparable. 

Our physical theory describes how helicase unwinding velocity and processivity depend on the interaction 
potential between the helicase and the ss-ds junction. This type of theory includes more detail then 
models based on kinetic states and transitions between them, while neglecting many details present in 
all-atom simulation models. Physical descriptions of the type we describe here complement other modeling 
frameworks. Simplifying some aspects of the system — for example, our choice to neglect the biochemical 
states which occur during helicase ss translocation — allows us to focus on the coupling between helicase 
translocation and NA unwinding. This level of detail may allow one to relate observed helicase behavior 
to simple physical mechanisms, and thereby gain more insight into the parameters which are important for 
helicase behavior. For example, experiments on UvrD suggest that this helicase unbinds more rapidly when 
unwinding dsDNA that during ss translocation Our model gives a physical picture for why this might 
happen: an interaction potential corresponding to active opening naturally leads to accelerated helicase 
unbinding. 

Our simple description neglects several effects. We ignore deformations of the NA strand, such as 
bending and torsion, and treat the strand as a rigid structure. The helicase is described by forward and 
backward rates only; we neglect the details of the protein's biochemical states. In addition, we ignore the 
effects of the NA base sequence on opening. These effects are believed to be weak for most helicases pp, 
although sequence dependence has been demonstrated for Rho 61 and DnaB 56 helicases. Recent work by 
Kafri, Lubensky and Nelson||)2] shows that a motor protein which translocates on a random track can show 
interesting behavior near the stall force. 
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